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Abstract 



We reanalyze the recent data from the Liquid Scintillator Neutrino Detector 
(LSND) experiment, that might indicate v'^ <-> v'e mixing. This indication is 
not completely excluded by the negative results of established accelerator and 
^^ , reactor neutrino oscillation searches. We quantify the region of compatibility 

QQ I by means of a thorough statistical analysis of all the available data, assuming 

0^ ■ both two-flavor and three-flavor neutrino oscillations. The implications for 

fvq ' various theoretical scenarios and for future oscillation searches are studied. 

^D ■ The relaxation of the LSND constraints under different assumptions in the 

statistical analysis is also investigated. 
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I. INTRODUCTION 

The Liquid Scintillator Neutrino Detector (LSND) |I| collaboration at the Los Alamos 
Meson Physics Facility (LAMPF) claim evidence for neutrino oscillations z/^ <-^ z/g from 
muon decay at rest 0,0], confirming preliminary results reported in 0]. 

The neutrino energy range, E,y ^ 20-53 MeV, and the path length, L ~ 30 m, make the 
LSND experiment particularly sensitive to neutrino mass differences in the eV range, which 
is of great interest for cosmological structure formation (see, e.g., 0). Mass differences of 
the same order are probed by other accelerator and reactor neutrino oscillation experiments 
^, which, however, do not find evidence for oscillations. Therefore, the interpretation of 
the LSND results as a neutrino oscillation effect requires careful checks and independent ex- 
perimental confirmation. In the meanwhile, it is legitimate to study its possible implications 
in the phenomenology of neutrino masses and mixing. 

The usual and easiest way to study the LSND signal is in terms of pure two-flavor 
oscillations. More precisely, the tau neutrino is assumed to be decoupled, so that only 
^fi ^^ '^e oscillation^ may occur. The LSND oscillation signal can then be compared, 
in the usual mass-mixing plane, to the limits coming from experiments probing the same 
flavor appearance channel (E776 and the Karlsruhe Rutherford Medium Energy Neutrino 
(KARMEN) experiment P|-pi]|]), as well as those probing z/g disappearance (reactors Gosgen 
|TT| , Bugey [|12|, and Krasnoyarsk |T^) and z/^ disappearance (CERN experiment CDHSW 
|T^). So far, this comparison has been simply performed by superposing the published 90% 
C.L. limits (see, e.g.. Fig. 31 in 0). However, the intersection of 90% C.L. regions is not 
a 90% C.L. region, and a more refined combination of the different experimental results is 
desirable. 

A three-flavor oscillation framework is required if, in addition, one wants to study the 

interplay between LSND and the tau neutrino appearance experiments, such as E531 at 

Fermilab [|l^ (completed), the CERN Hybrid Oscillation Research Apparatus (CHORUS) 

I6| , and the Neutrino Oscillation Magnetic Detector (NOMAD) experiment [0. The most 



general three-flavor approach would imply a reanalysis of all the data in a huge parameter 
space. However, if two of the three neutrinos are assumed to be almost degenerate in mass 
(so that the oscillations are driven by just one mass parameter at the laboratory scale), 
then the three-flavor analysis becomes much more manageable (see []T^,|TP[ and references 
therein). Many authors p0|-|30| have recently compared the LSND results with previous 



laboratory neutrino oscillation data within such approximation. In all these analyses pQ-pO 
the quoted limits on the neutrino masses and mixings were simply obtained as intersections of 
the region of parameter preferred by LSND and the regions excluded by the other oscillation 
experiments. As for the two-flavor case, no definite confidence level can be attached to such 
limits and to the corresponding scenarios, with the additional complication that the degrees 
of freedom increase when passing from a two-flavor to a three-flavor oscillation scheme. 
Motivated by the wide interest in the LSND results, and by the lack of a rigorous com- 



^In the following, neutrinos and antineutrinos will not be distinguished. Their oscillations in 
vacuum are equivalent both in the two-flavor limit and in the three-flavor scheme that we will use. 



parison with the estabhshed accelerator and reactor neutrino oscillation limits, we set out to 
quantify this comparison in a statistically consistent approach. We build upon our previous 
work |]18[, in which the results of the most constraining negative searches for neutrino oscilla- 



tions were reanalyzed ab initio and combined in a global x^ analysis. In the present work we 
include in the data set the recent LSND and KARMEN results, using both a two-fiavor and 
a three-flavor approach. We obtain precise bounds on neutrino mixings and mass differences, 
in the hypothesis that the LSND signal can indeed be interpreted as evidence for neutrino 



oscillations. These bounds are then used to discuss various scenarios described in pOj-pO 
as well as to investigate possible expectations for future laboratory oscillation searches. As 
a byproduct of our reanalysis of the LSND data, we gain interesting information on the 
stability of the neutrino oscillation fit. In particular, we show how the LSND constraints 
get weakened when the energy distribution of the LSND signal is gradually integrated out. 
Our work is organized as follows. In Sec. II we combine the LSND results with all the 
other established accelerator and reactor data relevant for two-fiavor oscillations z/^ ^-> u^, 
and derive bounds on the neutrino mass and mixing parameters. In Sec. Ill we include also 
the data from the established Ur appearance searches in a three-flavor oscillation framework 
with hierarchical mass differences. In Sec. IV we investigate the implications of our results on 
some phenomenological scenarios and on future oscillation searches. We draw our conclusions 
in Sec. V. Our reanalysis of the LSND data is presented in Appendix A, together with a 
critical investigation of the stability of the fit. The reanalysis of the KARMEN data is 
detailed in Appendix B. 



II. ANALYSIS IN TWO FLAVORS 

Let us assume pure two-fiavor oscillations between the electron and muon neutrino, with 
mass-mixing parameters (m^, sin^ 29). The event excess claimed by the LSND collaboration 
constrains these parameters in a relatively narrow region, as shown in Fig. 31 of P|. 

We have independently reanalyzed the published LSND data with a maximum likelihood 
method, as detailed in Appendix A. The results are shown in Fig. |l|. The 90 and 99 % C.L. 
contours compare well with the LSND confidence limits (Fig. 31 of 0), except for a "wiggle" 
at rri^ ~ 3 eV^ that, in our fit, extends to higher values of sin^ 29. This difference is irrelevant, 
being located in the region excluded by other established neutrino oscillation searches (as 
we shall see). Our best fit is reached for {m^, sin^ 26') = (3 eV^, 0.009). The best-fit point 
is not particularly important in itself, since good fits are also reached in other zones that 
form a set of almost-degenerate maxima of the likelihood function C, as can be realized by 
looking at the 68% C.L. contours of Fig. 1. One of the local maxima practically coincides 
with the best-fit point of the published LSND analysis (m^, sin^ 2^) = (19 eV^ 0.006). 
The fit disfavors values of m"^ close to integer multiples of 4.3 eV^, for which the oscillation 
probability of the high-energy neutrinos (which are responsible for most of the event excess) 
is suppressed. 

The reader is referred to Appendix A for an extensive discussion of the LSND data 
fit. Here we just anticipate that variations in the statistical treatment of the data may 
produce nonnegligible changes in the fit. In particular, the evidence in favor of neutrino 
oscillations appears to be weakened as the energy distribution of the LSND event excess 



is gradually integrated. Therefore, we think that the exact shape of the present LSND 
confidence contours should be taken with a grain of salt, while waiting for higher statistics 
and more stable fits. For definiteness, however, we will always refer to the LSND limits as 
derived from our maximum likelihood reanalysis (Fig. |l|) of the published LSND data. The 
combination with the data from other experiments is performed through the x^ method. To 
this purpose, the LSND likelihood is transformed as x^ = ~ In '^/2; the LSND x^ is then 
summed to the x^'s of the other experiments. 

The LSND indications can be directly compared with the (negative) results from two 
experiments that probed the same oscillation channel v^ <-> Ve-, namely KARMEN [^ and 
E776 0. Our reanalysis of the KARMEN data is described in detail in Appendix B. The 
results agree very well with the published oscillation bounds 0. Figure |^(a) shows our 
estimated 90 and 99 % C.L. contours for the KARMEN experiment (Ax^ = 4.61 and 9.21 
for two degrees of freedom) . 

Figure ^(b) shows the results of a combined fit to the LSND and KARMEN data. The 
regions allowed at 90 and 99 % C.L. are only slightly restricted with respect to Fig. |l], as most 
of the LSND oscillation signal is beyond the present KARMEN sensitivity. However, the 
KARMEN experiment should be able to test a significantly larger region of the parameter 
space in about two years 0. 

The results from the E776 experiment are more constraining than those from KAR- 
MEN. In Fig. ^(c) we report the E776 exclusion contours as derived by our reanalysis of the 



E776 data 1T||. The combination with the LSND and KARMEN data is shown in Fig. |(d). 
Notice that the 90% C.L. region of Fig. |^(d) is significantly larger than the intersection (not 
shown) of the 90% C.L. regions allowed by the LSND, KARMEN, and E776 experiments 
separately, especially for small mixing. The reason is that the E776 (and KARMEN) nega- 
tive results are compatible with no mixing, and thus push the combined fit towards values 
of the mixing angle smaller than those allowed by LSND alone. Moreover, the addition of 
the E776 data weakens the global evidence for oscillations, so that the no oscillation limit 
is allowed at 99% C.L. in Fig. 0(d), while it is excluded in Fig. ^(b). In other words, by 
adding the E776 and KARMEN data, the gradient of the x^ is suppressed for small values of 
the oscillation parameters, so that the 90% C.L. constraints are significantly relaxed in that 
direction, and the 99% C.L. left contour disappears. These results show that a combined 
analysis of the neutrino oscillation data may give rather different bounds than the naive "su- 
perposition" of oscillation limits, especially when the individual experimental results pull 
the fit in different directions. 

In the two-flavor approximation, the v^ disappearance searches at reactor experiments 



can be interpreted as probes of u^ ^^ ^a oscillations. The data from the Gosgen |TT[ , Bugey 



1^ , and Krasnoyarsk [|T3| experiments were reanalyzed and combined in [|T8|. The results 
are shown in Fig. ^(e). When combined with the LSND, KARMEN, and E776 data, the 
reactor results exclude the region at large mixing and low mass differences, as shown in 
Fig. 11(f). Notice that the narrow, disconnected regions allowed at 90% C.L. at large rn? in 
Fig. ^(f ) would disappear if a simple superposition of oscillation plots were made instead of 
a global fit. 

Finally, we mention that the results from v^ disappearance searches can also be added 
straightforwardly in the two-flavor approximation. However, the available data (CDHSW 
experiment []I1|) probe relatively large mixing angles, and we have checked that their addition 



produces negligible changes (not shown) to the contours of Fig. 0(f). The CDHSW data 



are anyway added, in the way described in ||T8[, in the global three-flavor analysis of the 
following section. 

III. ANALYSIS IN THREE FLAVORS 

In this section we combine with the x^ method the data from eight (reanalyzed) neutrino 
oscillation experiments (LSND, KARMEN, E776, Gosgen, Krasnoyarsk, Bugey, CDHSW, 
and E531 |]^), working in the simple three- flavor framework described in |T^. This frame- 



work is characterized by the assumption that two out of three neutrino mass eigenstates 
Vi {i = 1,2,3) are effectively degenerate in mass, Im^ — mfl ~ 0, while the mass gap of 
the "lone" state z/3 is taken in the range of laboratory neutrino oscillation experiments: 
m'^ = \m'^ — m\ ^| ^ 10~^ eV^. The neutrino oscillation parameter space is then spanned by 
three variables only: the dominant square mass difference m?, and the mixing angles and 
%l) (equal to ^13 and ^23 in the standard ordering of mixing matrices 0). More precisely: 

2|2 2|^^|2 2| /i\ 

m = \m^ — 1712 1 1 ^ l"^2 ~ ""^1 1 5 (l^J 

z/3 = sin0 z/e + COS0 (sin-i/; z/^ -|- cosz/' i^r) , (lb) 

with ip, (p & [0, vr/2]. We prefer to use the variables (tan^ ip, tan^ 0), which prove very useful 



for graphical representations ||18[. Notice that the two-flavor oscillation limits z/^ <-h> u. 



z/g <-!> Ur, and z/g ^^ z/^, are reached for = 0, ■?/' = 0, and ijj = 7r/2, respectively. The reader 



is referred to [IS] for further details and bibliography. 

Figure ^ shows the results of our analysis in the (?7i^, tan^ 0) plane, for twelve repre- 
sentative values of tan^ tp. The 90 and 99 % CL. contours shown in each panel represent 
sections (at fixed tan^ ^) of the confidence volume defined by x^ ~ Xmin = 6.25 and 11.36 
respectively (for three degrees of freedom). 

In the first panel of Fig. |^ the value of tan^ ip is very large, corresponding to almost 
pure u^ ^-^ Uf, oscillations {^p — * 7r/2). In fact, the limits close to the left side of the first 
panel represent a mapping of the two-flavor oscillation limits shown in Fig. |^(f) (modulo 
the different number of degrees of freedom). The mirror limits on the right (0 -^ 7r/2 — 0) 
in the same panel come from enlarging the mixing angle range from [0, 7r/4] (2 flavors) to 
[0, 7r/2] (3 flavors or more). 

In the subsequent panels the value of tan^ ip gradually decreases and the mixing of u^ with 
i>T- increases. The 2z/ left-right symmetry of the plots (0 — ^ 7r/2 — 0) is broken. The allowed 
regions get gradually shrinked, since for -0-^0 one approaches the two-flavor oscillation 
limit z/g ■^^ Vri which is not compatible with the LSND data. The allowed regions on the 
left half of each panel disappear more rapidly, since for t/) ^ and — > the relevant mass 
eigenstate (z/3) tends to become a flavor eigenstate (z/,-) and the oscillation phenomenon 
vanishes. Notice that the combination of all the data constrains m? above ~ 0.2 eV^ at 
90% CL. However, at 99% CL. there is no lower bound on m^ and the global fit becomes 
compatible with no oscillations. 

Figure ^ shows an alternative representation of the same confidence volume of Fig. |^. 
The volume is shown in (tan^'0, tan^ 0) sections at fixed values of m?. We remind that 



the left, right, and lower side of each panel correspond, asymptotically, to pure v^, ^^ z/,-, 
^e ^^ ^'/i; and z/^ <-s> Ut- oscillations, respectively. In Fig. H, the region allowed at 90% 
CL. generally consists of two disconnected parts at low and high 0. These zones would 
be continuously connected by a C-shaped band if only the LSND data were fitted 



However, the data from all other oscillation experiments strongly disfavor the central part 
of each panel, corresponding to large (and not observed) three-flavor mixing. The width of 
the 90% allowed regions is very sensitive to small variations of rri^ around multiples of 4.3 
eV^ where, as noticed, the LSND oscillation probability is suppressed. This is particularly 
evident in the panel at m^ = 13 eV^ where no combination of mixing angles is allowed at 
90% CL. As in Fig. |^, the no oscillation scenario (corresponding to the lower corners and 
the upper side of each panel) appears to be allowed at 99% CL. by the combination of all 
the data. 

Figures ^ and ^ are the main result of this work. They represent a concise summary (in 
a three-flavor scheme) of the most constraining neutrino oscillation data available, in the 
hypothesis that the LSND signal can indeed be interpreted as a signal of neutrino oscillations. 

IV. PHENOMENOLOGICAL IMPLICATIONS 

In this section we study the implication of the results shown in Figs. ^ and ^ for various 
theoretical scenarios and for short-baseline neutrino oscillation experiments. 



A. Implications for theoretical scenarios 

Apart from the LSND indications, there are (older) hints of neutrino oscillations from 
the solar neutrino problem and from the atmospheric neutrino anomaly. These additional 
data have been recently analyzed, within the same three-flavor framework adopted in this 
work, in [^] for solar neutrinos and in [^ for atmospheric neutrinos. 

The results obtained in this work and in [0,0 indicate that the LSND and atmospheric 
data are difficult to reconcile at 90% CL. if the subdominant neutrino mass difference is 
called to solve the solar neutrino problem [0]. In fact, from the LSND analysis (Figs. ^ and 



we derive m'^ ^ 0.22 eV^, while from the atmospheric neutrino analysis |^3| we derive 
the incompatible bound m'^ ^ 0.1 eV^. These data could be marginally reconciled only 
by dropping the information provided by the multi-GeV event sample of the Kamiokande 
experiment, as observed in |]2^ (see also |2^). If no data are excluded, however, one cannot 



achieve with three-flavor oscillations a good flt of the solar, atmospheric, and LSND data 
at the same time. The addition of a fourth, sterile neutrino might reconcile all the data 



without exclusions (see, e.g., [27]), but we think that more robust experimental checks must 
be performed before adopting such an extreme interpretation. 

Let us consider now the interplay between the LSND and the solar neutrino data. As 
shown in Fig. ^, at any given m'^ there are two regions allowed at 90% CL. at small and large 



0. However, only the small region survives the comparison with solar neutrino data [^ 
since at large the solar neutrino deflcit becomes almost energy- independent, contrary to 
the present experimental evidence |^^. The small 0, large ip allowed solution of Fig. ^ would 



indicate that the "lone" neutrino state u^ is dominantly a u^, thus excluding a simultaneous 
hierarchy of masses and mixings. 

Many semi-quantitative three-flavor analyses of laboratory neutrino data cornered a 
third, marginal "90% LSND allowed region" at small values of and ip p0| , |23| , ^ , |26| , |30 



corresponding to z/3 dominantly coupled with z/,-. This possibility is appealing since small 
mixings are generally regarded as more natural. However, our quantitative analysis shows 
that such region does not appear at 90% CL. for any value m^ and, therefore, it is unlikely 
that the z/^ is the dominant flavor component of u^. 

The threefold maximal mixing scenario (see, e.g., [§S|)' corresponding to 
{tair^ip, tan^ 0) = (1, 1/2), is excluded at more than 99% CL. for all values of m^ in the 
LSND sensitivity range, as can be seen in Fig. ^. The quasi-maximal mixing scenario pro- 
posed in |2^ , corresponding to ■?/' ~ 7r/4 and small 0, also appears to be strongly disfavored 



by our analysis. 

In conclusion, if one believes in the indications from solar and LSND neutrino experi- 
ments, then the dominant flavor component of the mass eigenstate u^ is the muon neutrino. 
This conclusion is relevant for model building. The precise bounds on the z/3 flavor content 
can be derived from Fig. ^ at any given value of m^, taking into account that the "large 0" 
solution is excluded by the analysis of solar neutrino data 



B. Implications for short-baseline experiments 

The implications of the LSND results for experiments probing the same oscillation chan- 
nel (z/^ <-^ z/e) are straightforward: they should observe a positive indication for neutrino 
oscillation, or disprove it, if they reach (at least) the LSND sensitivity. This goal should be 
reached by the KARMEN experiment in about two years 0. We mention in passing that 
the NOMAD experiment is also sensitive to Ue appearance for m'^ ^ 10 eV^, where it is 
expected to improve the existing hmits |^ . 



The implications for experiments probing different oscillation channels are less evident 
and require a three-flavor language. Here we focus on z/^ ^^ Uj. searches, currently being 
performed by the CHORUS [|1^] and NOMAD [^ experiments at CERN. Both experiments 
are expected to release soon the preliminary results of the flrst two-year run (1994-95), and 
are scheduled for additional two years of data taking (1996-97) |^. Two proposals for their 
upgrade are currently being considered at CERN, namely, the Tracking and Emulsion for 
Neutrino Oscillation Research (TENOR) [^, and the Neutrino Apparatus with Improved 
Capabilities (NAUSICAA) [^. The Cosmologically Signiflcant Mass Oscillation Search 
(COSMOS) experiment |^ at Fermilab is also planned to probe the z/^ <->■ z/,- channel 



with an expected sensitivity greater than CHORUS and NOMAD but somewhat lower than 
TENOR or NAUSICAA. 

Figure ^ shows the 90% CL. regions that can be probed by CHORUS or NOMAD (in 



two years |6|,|17l and four years i|,|3), COSMOS, and TENOR or NAUSICAA, using the 
{tan^ip, tan^ 0) representation at flxed values of m^. Also shown is the region preferred at 
90% CL. by the combination of all available data (including LSND), as taken from Fig. ^. 
It can be seen that the CHORUS and NOMAD experiments cannot probe, with a two-year 
statistics, the region preferred by all the data at 90% CL. (except for a marginal zone 



at m^ ~ 10 eV^), as already noticed in |T^. However, it is interesting to note that two 



additional years of data taking will make the CHORUS and NOMAD experiments able to 
probe a fraction of the small 0, large ip allowed region (see, e.g., the two top panels in Fig. ^. 
The COSMOS experiment and the TENOR or NAUSICAA projects will not only improve 
the sensitivity to the small solution at low m'^, but could even probe the large (p solution 
at large m^ (as in the first panel of Fig. ^. 

It follows from Fig. ^ that, if an oscillation signal shows up in CHORUS or NOMAD 
(after a 4-year run) and if the LSND claim is confirmed, then the mixing angles will be tightly 
constrained in the small 0, large ip region. Such a solution would be compatible with solar 



neutrinos ||3^, but not with atmospheric neutrinos I^Sj. Conversely, if the CHORUS and 
NOMAD oscillation searches give negative results, then the present small 0, large ip solution 
will be strongly reduced. The proposed COSMOS, TENOR, and NAUSICAA experiment 
might probe, in part, the additional possibility offered by the large solution. We conclude 
that the running and future u^ <-> 1/^ experiments can explore an interesting fraction of 
the LSND allowed region, and will provide in any case decisive constraints on the neutrino 
mixing angles. 



We finally mention that the Ug — >■ Ug reactor experiments in construction at Chooz 41 



and Palo Verde P2| are planned to improve the existing limits in the range 10~^ eV^ ^ 
m^ ^ 10~^ eV^. However, they are not expected to improve the available Bugey bounds in 
the range of m^ relevant for LSND (m^ ^ 0.22 eV^) and, therefore, should not have enough 
sensitivity to probe the 90% C.L. regions shown in Fig. ^. More precisely, the Chooz or Palo 
Verde sensitivity region would appear as a horizontal band in the (tan^ -0, tan^ 0) plane 



with an expected width given by 0.04 < tan^ ^ 25, which does not overlap with the 90% 
C.L. allowed regions of Fig. |[ 

V. CONCLUSIONS 

We have studied the interplay between the recent LSND indications and the negative re- 
sults of the established oscillation searches at accelerators (E776, KARMEN, E531, CDHSW) 
and reactors (Gosgen, Bugey, Krasnoyarsk). A thorough analysis has been performed as- 
suming two-flavor and three-flavor mixing. The region of neutrino oscillation parameters 
preferred by all the data has been determined. In particular, the three-flavor analysis in 
the (m^, ip 0) space typically favors two solutions, one at small and the other at large 0. 
Only the small solution is consistent with the solar neutrino data. A small fraction of 
the small solution appears to be in the range explorable by the CHORUS and NOMAD 
experiments (four year run). A larger fraction of this solution, and even a small part of 
the large region, could be probed by the COSMOS experiment at Fermilab and (with 
greater sensitivity) by the TENOR or NAUSICAA experiments at CERN. Finally, it has 
been shown (Appendix A) how different statistical treatments of the published LSND data 
affect the bounds on the neutrino oscillation parameters. 
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APPENDIX A: REANALYSIS OF THE LSND DATA 

The Liquid Scintillator Neutrino Detector (LSND) experiment [|l| searches for z/g appear- 
ance in a Uf^ beam from the decay (at rest) 

-^ e+ + i^e + «>/. , (Al) 



through the reaction 



Ue+P ^ e^ +n . (A2) 



The maximum z/^ energy is E^'"^^ = 52.8 MeV. For any candidate positron, the LSND 
detector can determine both its energy Ee and its vertex coordinate L (equal to the neutrino 
path length). The detector covers the range L = 29.8 ± 3.8 m. 

Positron-like events have been observed 0,^ in excess of the estimated background. If 
this excess is interpreted as a signal of neutrino oscillations, a region of (two-flavor) neutrino 
mass and mixing parameters appears to be preferred by the data analysis, as shown in 
Fig. 31 of i. 

In this Appendix we describe how such region can be reproduced to a good accuracy by 
reanalyzing the published LSND data, and how it gets modified by changing some assump- 
tions. We fix the terminology in Sec. Al, develop the standard analysis in Sec. A2, and 
discuss variations of the standard analysis in Sec. A3. 

1. Basic definitions 

The total LSND signal 5* is defined as the sum of the total (beam-off plus beam-related) 
background B and of the event excess E due to possible oscillations, 

S = B + E . (A3) 

The subscript "exp" and "theo" will be used to distinguish the observed values of S 
and E from the theoretical estimates, respectively. No subscript is attached to B, which is 
always a simulated quantity. 



For a given experimental distribution of the signal (divided in N bins), {ni}i=i„,N, and 
the associated theoretical distribution {fii}i=i,,,N, a likelihood function £ can be defined as 



N 



'^ = n^/^*"'e-'^" • (A4) 



• 1 n, 



The analysis performed by the LSND collaboration is essentially based on the likeli- 
hood Cel associated to the (experimental and theoretical) double differential distribution 
(PS/dEe dL, 



-EL 



C{d^S/dE^dL) , (A5) 



having divided the positron energy range (20 MeV < -Eg < 60 MeV) into 20 bins and the 
length range (26.0 m < L < 33.6 m) into 38 bins.0 

Unfortunately, the double differential distributions d^Sx/dEe dL {X = exp, theo) are not 
published in Refs. [00, so that the likelihood function in Eq. ( [A5| ) cannot be reproduced 
exactly. However, one can recover at least the projected distributions dSx/dEe and dSx/dL 
(X = exp, theo) from [^,0 and thus calculate the corresponding likelihoods Ce and Cl- In 
the absence of more detailed information, we assume that the total likelihood is factorizable, 

Cel ^ Ce X jCl 

= C{dS/dE^) X C{dS/dL) , (A6) 

as if dS/dEe and dS/dL were independent. Although such statistical independence is never 
exactly realized, we will see a posteriori that the published LSND oscillation analysis 0] is 
reproduced rather well through the approximation Eq. ([A6| ) . 

In the next subsection we discuss how the relevant distributions of the signal can be 
recovered from Refs. [|l],0 and then used to perform a maximum likelihood analysis. 

2. Standard Analysis 



The approximate likelihood function in Eq. ( |A^ ) is based on the four distributions 
dSexp/dEe, dSexp/dL, dStheo/dEe, and dStheo/dL. The first two can be read off Fig. 30(a) 
and 30(d) of respectively (dots with statistical error bars). The last two are estimated 
as follows. 

The double differential (theoretical) distribution of events E in excess of the background 
is calculated as 

d^Ethco .r /■ ,p, d^^ 2 



dE.dL 



M jdK -^ L-' <E^) <K) R{Ee, K) PiL/E^) , (A7) 



^Actually, the complete LSND analysis H includes also subdominant information related to the 
positron scattering angle and the event pattern recognition. 
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where E'^ and Ee are the true and measured positron energy respectively, E^, is the neutrino 
energy, d^^/dE^ is the neutrino energy spectrum, a is the neutrino cross section, e is the 
detector efficiency, R is the resolution function, P is the oscillation probability, and TV is a 
normalization factor. 

The Michel energy spectrum d^v/dEy of z/^'s from /x^ decay at rest is well known (see, 
e.g.. Fig. 6 in [|l|). The total flux varies approximately with the inverse square of the distance 
L (see Sec. II B of 0^ 



The neutrino cross section a for the reaction ([A2| ) is also well known [^. The neutrino 
energy and the (true) positron energy E'^ are tightly related, E^ = E'^ + 1.8 MeV + 5{6e)i 
where 5 is a small (and often neglected) kinematical correction depending on the positron 
scattering angle Oe- We have applied a fixed correction corresponding to the average scat- 
tering angle cos6'e ^ 0.2 (see Sec. VI C of 0). 

The efficiency s{E'^ for detecting a positron with true energy E'^ can be recovered from 
Fig. 9 in [g]. It varies from 40% to 90% above the analysis threshold {E^ > 20 MeV). 

The difference between the measured and true positron energy is well approximated by a 
Gaussian energy resolution function R{Ee, E'J with one-sigma width equal to 7.7% at 52.8 
MeV, and scaling as l/E'^ for other energies (see Sec. VI D of IQ). 

Finally, the normalization factor Af is fixed by imposing that for 100% transmutation 
(P = 1) the event excess -Etheo (integrated over the detector length and over the energy 
range 20 MeV < E^, < 60 MeV) be equal to 16670 events, as reported in Table IV of 0. 

We have checked that, for P = 1, the shape of dEtheo/dEe derived from Eq. ([A7| ) compares 
very well with the corresponding LSND simulation (Fig. 7 in 0). Moreover, our estimated 
value of -Etheo, integrated over the energy subrange 36 MeV < Ee < 60 MeV, coincides with 
the LSND quoted value (12500 events, see Table III of 0) within 0.5%. These two nontrivial 
checks add confidence in our calculation of the event excess due to neutrino oscillations 
through Eq. {\K7\ }. 

Concerning the background B, the relevant distributions dB/dE^ and dB/dL are not 
reported explicitly in the LSND publications [|l],0]. We have derived them indirectly by 
subtraction, B = S — E. More precisely. Figs. 30(a) and 30(d) in p show (as solid lines) 
the distributions dStheo/dE and dSthco/dL of the signal expected for the oscillation case 
{Ami, sin^'^d) = (19 eV^ 0.006). We calculate dEtheo/dE and dE^^^^/dL from Eq. (^) for 
the same oscillation parameters and then obtain the background distributions by subtraction, 

dB = dStheo — dEtheo- 

Figs. ^ and show, respectively, the distributions dB/dE^ and dB/dL derived in this 
way (dashed histograms). The integrated background is P = 1701 events. To guide the eye, 
we also show in Figs. ^ and |7] the corresponding distributions of Sexp (dots with error bars) 
and of 5theo [for (m^, sin2 2e) = (19 eV^ 0.006)], as taken from Figs. 30(a) and 30(d) of 
respectively. The integrated signal is ^exp = 1763 events, corresponding to S'exp — P = 62 
excess events. 

With the above ingredients we calculate, for a given oscillation probability function P, the 
E and L distributions of the expected signal ^theo and the associated likelihood [Eq. (^)], 
from which the confidence contours of Fig. |l| have been obtained. 

When used in combination with other neutrino oscillation data, the LSND likelihood is 
transformed in a x^ through the relation f^ X^ = — ln£/2. 
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3. Variations in the analysis 

The "standard" analysis of the LSND data described in the previous section (and used 
throughout this work) makes use of the E^ and L distributions of the data (Figs. || and ^ 
through a maximum hkehhood analysis. 

In order to test the stability of the "standard" LSND limits shown in Fig. |l|, we have 
performed a few alternative analyses of the LSND data that gradually deviate from the 
standard one. The results are shown in Fig. ^ 

Fig. |](a) shows the 90 and 99 % C.L. limits obtained by dropping the likehhood factor 
Cl in Eq. (|A^). The limits are only slightly relaxed with respect to Fig. |T]. In fact, the fit 
is basically driven by the excess of events in a few bins of the energy distribution (Fig. ^, 
while the path length distribution (Fig. |^) is not really discriminating within the statistical 
error bars. For the sake of simplicity, only the dominant E^ distribution is used in all panels 
of Fig. §. 

Fig. ^(b) shows the LSND limits from a "true" x^ analysis of the E^. distribution (i.e., 
the canonical definition of x^ is used and not — lnC/2). Only statistical errors are included. 
A comparison with Fig. ^(a) shows a slight relaxation of the C.L. contours for small mixing 
(i.e., for lower rates), as expected from the Gaussian approximation of a Poisson distribution 
which is implied by the x^ method. 

Fig. H(c) shows the results of the same x^ analysis of Fig. ^(b) with the addition of a 
plausible 10% uncertainty in the background normalization 0. The correlation of this error 
between any two bins is taken equal to 1. The LSND oscillation limits are somewhat relaxed, 
and the data became even compatible with no oscillations at about 99% C.L. 

In Fig. H(d) the same analysis of Fig. |](c) is repeated by dividing the energy distribution 
in 5 bins instead of 20. The rationale for this exercise is the erratic position of bins where 
there is a significant event excess (see Fig. ^, that might be a symptom of statistical 
fluctuations. Such fluctuations should be somewhat flattened by grouping bins. In fact, 
the LSND limits appear to be slightly relaxed as compared with Fig. ^(c). A more evident 
effect is obtained by dividing the Eg distribution in just two bins, as shown in Fig ^(e). 

Finally, Fig. ^(f) shows the LSND limits obtained when the E^ distribution is fully 
integrated, i.e., when only the total number of events is used. Since the total event excess 
(62 events) is smaller than the systematic uncertainty of the background (10% of 1701 
events), no significant indication in favor of neutrino oscillations is obtained, and the analysis 
gives only exclusion contours. Therefore, the preference for nonzero values of the neutrino 
oscillation parameters in Figs.H(a)-(e) appears to be driven by the information contained in 
the positron energy spectrum. 

In conclusion, the exercises described in this section show that, as far as the published 
information is used, the "standard" LSND limits of Fig. |I] appear to be dominated by the 
detailed energy distribution of the observed event excess. Different ways of treating this 
distribution may lead to significant changes in the C.L. contours. In particular, the indica- 
tions in favor of neutrino oscillations are increasingly weakened by gradually integrating the 
energy spectrum information. 
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APPENDIX B: REANALYSIS OF THE KARMEN DATA 

The Karlsruhe Rutherford Medium Energy Neutrino (KARMEN) experiment P|-pi]| is 
being performed at the pulsed spallation neutron facility ISIS of the Rutherford Appleton 
Laboratory. A search is made for z/g and v^^ appearance from i/^ and z/^ produced in 7r+ 
and /i+ decay at rest. Therefore, the energy spectrum of the neutrino source, d^^/dE, is 
the same as in the LSND experiment. The KARMEN detector, a 56 ton liquid scintillation 
calorimeter, is located at a distance L = 17.5± 1.75 m (front-end range) from the source. 
The pulsed time structure of the beam can be exploited to achieve a strong suppression of 
the background. 

At present, the number of candidate events is consistent with the estimated background 
for both oscillation channels, z/^ <-»• Ue and z/^ <-> Ug. These negative results can be used 
to constrain the neutrino mass and mixing parameters. Our statistical analysis of the con- 
straints for the two oscillation channels and their combination is detailed below. 



1. Channel f^ <-> i/g 

For this channel we use three basic inputs 0: (1) the number of events (prompt 
positrons) observed in excess of the expected background, A^cxp ± cr^v = —^-^^-f^] (2) the 
number of events expected for 100% transmutation, A^max = 3038; and (3) the positron en- 
ergy spectrum S{Ee+) (Fig. 6 of ||TU|) expected for 100% transmutation. The above inputs 
refer to the default energy window 20 < Ef,+ < 52 MeV. The approximate energy relation 
Ep^ ~ Ee+ + 1.8 MeV is adopted. 

The expected (theoretical) number of positrons for a generic oscillation probability P = 
P{L/EpJ is then given by 

rL+D/2 



/ dLL-^ dE,+ SiE,+ )P{L/E,^ 

JL-D/2 J 

rL+D/2 } 

/ dLL-^ dE^+S(E^+) 

Jl-D/2 J 



T\j AT J L—D/2 J ('R1 "\ 

-'''theo -'''inax /•J + D/2 c ' V/ 




where D is the detector length (3.5 m). 

The probability e that a difference as large as A = (Atheo — A'exp)/o"Af be a statistical 
fluctuation is obtained by integrating a one-sided Gaussian (Atheo being nonnegative), 

2 r^ 

dx exp(-xV2) . (B2) 

As a check, we draw in Fig. ^ the exclusion curve (dotted line) corresponding to 1 — e = 
90% in the usual two-flavor mass-mixing plane. This curve matches well the 90% CL. 
contour reported by the KARMEN collaboration in |^ for the same oscillation channel. 

When a combination with other oscillation data is performed in a global x^ analysis, we 
use a "fake" x^ statistic determined by probability inversion [x^ = x^(e)]; yielding for any 
assigned C.L. the same exclusion curves. 
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2. Channel u^ <-> fe 

In this channel the energy spectrum is monochromatic [Ey = 29.8 MeV). We use three 
basic inputs 0: (1) the number of candidate events, N = 1; (2) the estimated background 
A''^ and its 1-sigma uncertainty, A''^ ± o-^ = 2 ± 0.3; and (3) the number of events expected 
for 100% transmutation, A^max = 154. 

For a generic oscillation probability, the expected number of events is then calculated as 



f.L+D/2 

I 
I 

7L+D/2 



/ dLL-^P{L/E. 

JL-D/2 

dLL-^ 



-''theo -''max pT + Dl2 ' \ / 



L-D/2 

The confidence level 1 — e associated to A^theo can be calculated by using the statistics 
appropriate to a Poisson process with background 0, with an allowance for statistical 
fluctuations of the background itself p4| : 



^ 



driB fins) e-("-+^'-°) 5:(nB + N.^^XM 



l-e = l ^ , (B4) 

Jo Z^n 



where 



fi^B) = -j= — exp 



2na 



B 



n=0 



1 fU B-Ns ^ '^ 

ctb 



(B5) 



The curve corresponding to 1 — e = 90% in the usual mass-mixing plane is shown in 
Fig. H as a thin, solid line. It reproduces in detail the 90% C.L. contour reported by the 
KARMEN collaboration in §]. 

As for the antineutrino channel, we use a "fake" x^ statistic determined by probability 
inversion [x^ = X^i^)] when a combination with other data is performed. In particular, we 
combine the two KARMEN oscillation channels, z/^ ^^ Ue and z/^ ^^ Ue, by summing the 
corresponding x^'s- The curve corresponding to a variation of 4.61 in the total x^ (90% 
C.L. for two degrees of freedom) is shown in Fig. |^ as a thick, solid curve (the same curve 
reported in the top left panel of Fig. ^. It can be seen from a comparison of the three curves 
in Fig. ^ that the channel 9^ ^^ Ug dominates the KARMEN exclusion region. 
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FIGURES 

FIG. 1. Region of the 2u oscillation parameters preferred by the LSND data (our reanalysis) 
at 68, 90, and 99 % C.L. (TVdf = 2). 

FIG. 2. Laboratory oscillation data in two flavors (our reanalysis). Left panels: contours of 
regions excluded at 90 (solid) and 99 % (dashed) C.L. by the accelerator experiments KARMEN 
and E776, and by reactor experiments (Gosgen, Bugey, and Krasnoyarsk combined). Right panels: 
variations in the region preferred by the LSND data (see Fig. |^) with the progressive addition of 
the KARMEN, E776, and reactor data. 

FIG. 3. Three-flavor analysis of the most constraining laboratory oscillation experiments 
(LSND, KARMEN, E776, E531, CDHSW, Gosgen, Bugey, and Krasnoyarsk combined). The pre- 
ferred region in the 2>v parameter space (m^, tan^ ^, tan^ (p) is shown through twelve (m^, tan^ (/>) 
sections at fixed, representative values of tan^V- Solid lines: 90% C.L. contours (Ax^ = 6.25 for 
iVoF = 3)- Dotted lines: 99% C.L. contours (Ax^ = 11.36). 



FIG. 4. As in Fig. y, but in (tan^ (j), tan^ ■0) sections at twelve representative values of 



2 

m . 



FIG. 5. Regions of the parameter space explorable at 90% C.L. by the following t/^ — >■ Ut 
experiments (in order of increasing sensitivity): CHORUS or NOMAD in two years (thick, dotted 
line) and four years (thin, sohd line), COSMOS (dashed line), and TENOR or NAUSICAA (thin, 
dotted line). These experiments can probe a fraction of the zone preferred at 90% C.L. by the 
combination of all the available data, including LSND (thick, solid line). 

FIG. 6. Energy distribution dS/dEe of the LSND signal (20 bins). Dashed line: 

background component dB/dEg (our reanalysis). Solid line: signal dSthco/dEe expected for 
{m?, sin^ 20) = (19 eV , 0.006), as taken from Fig. 30(a) of Q. Dots with statistical error bars: 
observed signal dSexp/d-Eej from Fig. 30(a) of 0. 

FIG. 7. As in Fig. |, but for the path length distribution dS/dL of the LSND signal (38 bins). 

FIG. 8. Variations in the LSND bounds with respect to the "standard" bounds of Fig. Q, as a 
result of of data analyses different from Cel maximization, (a) Maximization of Ce only, (b) x^ 
analysis of the energy distribution of the signal with statistical errors only, (c) x^ analysis of the 
energy distribution, assuming a systematic 10% uncertainty in the overall background normaliza- 
tion, (d) As in (c), but dividing the energy distribution in 5 bins, (e) As in (c), but dividing the 
energy distribution in 2 bins, (f) As in (c), but integrating the total signal (= 1 bin). See the text 
for details. 

FIG. 9. Results of our reanalysis of the KARMEN data for the neutrino and antineutrino 
channels and their combination. Contours are drawn at 90% C.L. (A^^ = 4.61 for A'^df = 2). 
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2v oscillations 




FIG. 1. Region of the Iv oscillation parameters preferred by the LSND data (our reanalysis) at 
68, 90, and 99 % C.L. (A^df = 2). 
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FIG. 2. Laboratory oscillation data in two flavors (our reanalysis). Left panels: contours of 
regions excluded at 90 (solid) and 99 % (dashed) C.L. by the accelerator experiments KARMEN 
and E776, and by reactor experiments (Gosgen, Bugey, and Krasnoyarsk combined). Right panels: 
variations in the region preferred by the LSND data (see Fig. |^) with the progressive addition of 
the KARMEN, E776, and reactor data. 
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3]y oscillations 

LSND + KARMEN + E776 + E531 + CDHSW + Goesgen + Bugey + Krasnoyarsk 
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FIG. 3. Three-flavor analysis of the most constraining laboratory oscillation experiments 
(LSND, KARMEN, E776, E531, CDHSW, Gosgen, Bugey, and Krasnoyarsk combined). The pre- 
ferred region in the 3u parameter space (m^, tan^ ^, tan^ cj)) is shown through twelve {m?, tan^ (/)) 
sections at fixed, representative values of tan^ tp. Solid lines: 90% C.L. contours (Ax^ = 6.25 for 
iVoF = 3). Dotted lines: 99% C.L. contours (Ax^ = 11.36). 
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3iy oscillations 

LSND + KARMEN + E776 + E531 + CDHSW + Goesgen + Bugey + Krasnoyarsk 
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FIG. 4. As in Fig. ^, but in (tan^ (p, tan^ -0) sections at twelve representative values of m^ 
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LSND + KARMEN + E776 + E531 + CDHSW + Goesgen + Bugey + Krasnoyarsk 

CHORUS/NOMAD, 2 yr (1994-95) COSMOS at Fermilab 

CHORUS/NOMAD, 4 yr (1994-97) TENOR/NAUSICAA at CERN 
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FIG. 5. Regions of the parameter space explorable at 90% C.L. by the following v^ -^ Vr 
experiments (in order of increasing sensitivity): CHORUS or NOMAD in two years (thick, dotted 
line) and four years (thin, sohd line), COSMOS (dashed line), and TENOR or NAUSICAA (thin, 
dotted line). These experiments can probe a fraction of the zone preferred at 90% C.L. by the 
combination of all the available data, including LSND (thick, solid line). 
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FIG. 6. Energy distribution dS/dEe of the LSND signal (20 bins). Dashed hne: background 
component dB/dE^ (our reanalysis). Sohd Hne: signal dStheo/dE^ expected for {m?, siv?26) = 
(19 eV , 0.006), as taken from Fig. 30(a) of [^. Dots with statistical error bars: observed signal 
dS^^p/dEe, from Fig. 30(a) of §. 
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dS/dL distribution 
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FIG. 7. As in Fig. ^, but for the path length distribution dS/dL of the LSND signal (38 bins). 
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Different 2v analyses of the E^ distribution (LSND) 
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FIG. 8. Variations in the LSND bounds with respect to the "standard" bounds of Fig. ||, as a 
result of of data analyses different from Cel maximization, (a) Maximization of Ce only, (b) y^ 
analysis of the energy distribution of the signal with statistical errors only, (c) x^ analysis of the 
energy distribution, assuming a systematic 10% uncertainty in the overall background normaliza- 
tion, (d) As in (c), but dividing the energy distribution in 5 bins, (e) As in (c), but dividing the 
energy distribution in 2 bins, (f) As in (c), but integrating the total signal (= 1 bin). See the text 
for details. 
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FIG. 9. Results of our reanalysis of the KARMEN data for the neutrino and antineutrino 
channels and their combination. Contours are drawn at 90% C.L. (Ax^ = 4.61 for Ndy = 2). 
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